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Molecular dynamics simulations in three-dimensional copper have been performed to quantify the 
void coalescence process leading to fracture. The correlated growth of the voids during their linking 
is investigated both in terms of the onset of coalescence and the ensuing dynamical interactions 
through the rate of reduction of the distance between the voids and the directional growth of the 
voids. The critical inter-void ligament distance marking the onset of coalescence is shown to be 
approximately one void radius in both measures. 
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The point at which voids begin to coalesce during dy- 
namic fracture is of considerable interest because com- 
plete fracture of the material typically ensues rapidly 
thereafter. A robust model of the onset of coalescence 
is an important ingredient of a predictive model of dy- 
namic fracture. The conventional picture of how ductile 
metals break under the rapid application of stress con- 
sists of three stages: void nucleation, growth and coales- 
cence. Initially voids nucleate from the weak points in the 
material such as inclusions and/or grain boundary junc- 
tions. Once nucleated, the voids grow under the tensile 
stress, driven by the reduction in elastic energy. Eventu- 
ally, the voids grow sufficiently large that they interact 
with each other, coalesce into larger voids, and finally 
form the fracture surface There are many interest- 

ing twists and subtleties, such as the interplay between 
shear localization and void growth, but the basic picture 
applies to the fracture of a broad class of ductile metals. 
Naturally considerable effort has gone into the study of 
this fracture process, both in modeling and theory and 
in experiment, including a new generation of 3D, non- 
destructive fracture characterization techniques such as 
x-ray tomography. Nevertheless, a robust, mechanistic 
understanding of coalescence has yet to emerge. 

Computationally void growth has been studied exten- 
sively at the continuum leveL ELS IE an d more re- 
cently at the atomistic level p], Isl 111 • The atomistic stud- 
ies demonstrate that voids grow by emitting dislocations, 
which carry away the material, platelets of atoms, from 
the void and are responsible for the plastic deformations 
needed to accommodate significant void growth. There 
are also many recent studies of fracture in ductile met- 
als with several holes or voids While these studies 
model the void growth explicitly, often with fairly so- 
phisticated models of plasticity, they typically simplify 
the coalescence process to instantaneous unification of 
the voids based on a relatively simple criterion such as 
growth of the voids to within one void diameter of each 
other or a plastic strain threshold. The earlier continuum 
studies (cf. Ref. and the one atomistic study known 
to us 11] of the coalescence process have been conducted 
in effectively two dimensional and highly symmetric sys- 



tems. 

In this Letter we analyze the details of the onset of void 
coalescence. In particular we quantify the point at which 
coalescence begins, as measured by a critical inter-void 
ligament distance (ILD), and examine the mechanisms 
involved in the transition from independent void growth 
to coalescence. There are several ways in which two voids 
can interact. In the case of pure impingement, the voids 
only interact when they grow to the point that they in- 
tersect and join into a single void. In reality, the voids 
interact before they intersect. Their range of interac- 
tion is extended due to their elastic and plastic fields. 
Each void generates an elastic strain field of the form 
generally associated with centers of dilatation 12]. The 
shear stress decreases with the distance from the void 
like r -3 . For voids sufficiently close each void's growth 
rate is altered by the stress field of the proximal void. 
The modification of the elastic field can affect the initia- 
tion of plasticity, as well as the subsequent development 
of the plastic zone around the voids. The voids may in- 
teract through their plastic fields, too, in which case the 
fields may give rise to an increased hardening rate in a 
localized region or to thermal softening and shear local- 
ization. An argument due to Brown and Embury for a 
transition to shear deformation based on simple geomet- 
rical considerations suggests that the critical inter-void 
ligament distance, ILD C , should be equal to one diam- 
eter of a void 0; that is, when the surfaces of a pair 
of voids are separated by one void diameter, they tran- 
sition from independent void growth to coalescence. It 
is at this point, they argue, that the dominant void pro- 
cess switches from the radial plastic flow around isolated 
growing voids to a shear deformation allowing the rapid 
coalescence of the pair of voids. However, more recent 
two-dimensional studies suggest that for distances be- 
tween voids as large as six diameters the void growth 
rate is enhanced 14]. 

The use of atomistic techniques permits the analysis of 
the contributions of these competing mechanisms to the 
onset of void coalescence, as we describe in this Letter. 
We demonstrate the existence of, and compute, the crit- 
ical inter-void ligament distance ILD C by starting with 
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two voids well separated from each other and detecting 
the point at which correlated growth begins, marked both 
by the accelerated rate at which the two void surfaces ap- 
proach each other and by biased growth causing the voids 
to start to extend toward each other. This gives an indi- 
cation of the onset of the coalescence process, and it tests 
the argument by Brown and Embury We also test 
the setup by Horstemeyer et al 14] by varying the initial 
distances between the voids and measuring the asymp- 
totic growth rate of the voids. The initial void-to-void 
distance below which the growth-rate is enhanced should 
give another candidate for the critical distance and mea- 
sure it in a volumetric sense. 

We have performed large-scale (parallel) classical 
molecular dynamics (MD) simulations in single crys- 
tal face-centered cubic (FCC) systems using an empirical 
embedded-atom model (EAM) potential for copper 
The three dimensional (3D) simulation box consists of 
120 x 120 x 120 4-atom FCC unit cells with periodic 
boundary conditions for a total of 6 912 000 atoms. The 
system is initially equilibrated using a thermostat at 
T = 300 K and a constant volume L 3 (with L = 43.4 nm) 
chosen to give ambient pressure, P c± MPa. Once the 
system has reached equilibrium, two spherical voids are 
cut in the system with radius ro = 0.05L = 2.2 nm: one 
in the middle of the box and the other 12.2 nm away, in 
a relative position of [0.25, 0.1166, 0.0544]L. We refer to 
these as void A and void B, respectively, see Fig. [I] When 
the initial distance between the voids is varied, the loca- 
tion of the void B is changed, but the relative orientation 
of the voids is kept fixed. Initially, the voids are equal in 
size, with approximately 3620 atoms removed for each. 
Once the voids are formed, the thermostat is turned off, 
and dilatational strain is applied uniformly at a constant 
strain-rate e. Applied strain-rates of e = 10 8 /sec and 
10 9 /sec have been used with perfectly triaxial, or hy- 
drostatic, expansion. The dilatation is simulated by ex- 
panding the simulation box at each time step 7], im- 
plemented in the fashion of the Parrinello and Rahman 
technique [HI in order to prevent the spurious generation 
of elastic waves at the box boundaries. More details of 
the simulation method can be found from Ref. 0. 

While some void growth takes place through elastic 
stretching in the initial phases of the box expansion, sig- 
nificant void growth and void- void interaction take place 
only once plastic deformation has begun. The impor- 
tant role of plasticity forces us to consider in some detail 
how dislocations are generated and the effect the disloca- 
tion dynamics have on void coalescence. Figure ^shows 
a visualization within a slice of width 4.5 A of a plane 
including centers of both voids at six different instants 
during coalescence. The atoms shown are either on the 
surface of the voids or belong to dislocations. The de- 
cision of which atoms to plot is based on a geometrical 
criterion, a finite-temperature generalization of the cen- 
trosymmetry deviation || . From the snapshots one 
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FIG. 1: Snapshots of the slices of the two- void system at e = 
10 9 /sec with only those atoms shown that are in dislocation 
cores, stacking faults, void surfaces or other defects (see text). 
The dashed loop in panel (c) is drawn around a slice of a 
prismatic dislocation loop. The plane shown passes through 
the centers of both voids. The snapshots show the initial 
plasticity (a,b), interacting plastic zones (c,d) and the final 
coalescence (e,f). The frames correspond to strains of e = 
1.72%, 2.42%, 3.47%, 3.89%, 4.52%, and 5.21%, respectively. 



sees that the deformation mechanism involves the nucle- 
ation and propagation of dislocations, accommodating 
the void growth, and the interaction of the dislocations. 
For example, the prismatic dislocation loops punched out 
by the voids appear as roughly parallel line traces (due 
to the stacking fault ribbons) in the slice Fig. [He), as 
verified in the full 3D configuration. Initially the dislo- 
cation activity around each void is essentially symmet- 
ric [Fig. [Ha) and (b)], as expected for independent void 
growth, but as the plastic fields evolve the void-void in- 
teraction is clearly evident both through interactions be- 
tween the two plastic zones and bias due to the elastic 
fields [Fig.njc)]. Once the dislocation density grows suf- 
ficiently high in the ligament region between the voids 
[Figs.[Hd)], v °id B begins to grow in the direction away 
from void A. Eventually the voids coalesce [Fig. [He)], 
and continue to grow as one until ultimately coalescence 
with the void's periodic images takes place [subsequent 
to Fig. Elf)] so that the cavity percolates through the 
periodic system. 

Figure |I] offers several visual indications of the inter- 
action between voids. Clearly, the separation between 
the void surfaces (the ILD) serves as something akin to 
a reaction coordinate for the coalescence: the voids coa- 
lesce when it goes to zero. Other indications include the 
displacement of the center of a void as it grows preferen- 
tially toward the neighboring void and the change in the 
void growth rate as the voids interact. We now quantify 
these effects in order to analyze the coalescence. 

In Fig. [3 the dynamic evolution of the ILD has been 
plotted for strain-rates e = 10 8 /sec and 10 9 /sec and 
for various initial closest surface-to-surface distances be- 
tween the voids ILDq. In Fig. El the case ILDq ~ 1.8 
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FIG. 2: (color) The dynamical ILD, the distance between 
the surfaces of the voids along the line connecting the orig- 
inal center positions, plotted versus the mean void size f 1 ^ 3 
for various initial ILDo's. The ILD is plotted in units of the 
average void diameter d calculated at each instant using the 
formula d = 2[3/(47r)/V] 1//3 assuming roughly spherical voids. 
For reference, thin lines are plotted to show the relationship 
for spherical voids impinging freely on each other (see text). 
The red line shows the hypothetical ILD computed by dupli- 
cating the single void at fixed centers with ILDo = 1.8. The 
horizontal line is at ILD=0.5d. 



was shown. While the strain is the quantity controlled in 
these simulations, we have found that the data collapse 
well and the coalescence is indicated much more clearly 
if the ILD data are plotted as a function of the linear 
mean void size, J 1 / 3 , where / is the average of the two 
void fractions / = V VO id/V and V is the volume of the 
box. The technique for calculating the void volume is 
described in Ref. 9]. The reason is that for voids grow- 
ing independently, the ILD is closely related to the size 
of the void through geometry. In particular, the ILD for 
two spherical voids of the same diameter d growing inde- 
pendently is given by ILD = (ILDo + do)(l + £) — d : where 
do is the initial diameter. This formula is used to gener- 
ate the thin curves in Fig. Deviations from these free 
impingement curves in the MD simulations indicate void 
shape changes, most importantly anisotropic growth due 
to coalescence. The same data plotted versus strain ex- 
hibit a strong strain-rate dependence inherited from the 
void growth law V VO id(s) (cf. Ref. |£j). For reference, the 
snapshots in Fig.^a)-(d) correspond to mean linear void 
size J 1 / 3 = 0.089, 0.094, 0.150, and 0.195, respectively. 
After coalescence / is not calculated. 

Initially the separation distance decreases essentially 
smoothly until the plasticity begins, eventually reaching 
zero. A transition occurs when the ILD starts to de- 
crease much faster than the free impingement line, which 
takes place when the ILD reaches approximately one half: 
ILD C = 0.5 ±0.1 diameter or one radius, independently 
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FIG. 3: Distance d cv from the original center of void to the 
instantaneous void center, projected onto the line connecting 
the original void centers, plotted versus the average void size 
to the point of coalescence (ILD ~ 0) . The sign of the distance 
d cv is positive for motion toward the other void. ILDo =1.8. 
The thin solid line is for a single void in the same size of the 
box and with the same radius and e — 10 9 /sec projected to 
the same line. Here the distance d cv is given in the units of 
the original void diameter do. 



of ILDo or the strain-rate. Note that the unit of ILD is 
the current diameter of a void, not the initial value. 
A curve derived from a single void growth is provided to 
estimate the contribution of uncorrelated faceting effects 
(the "one void" curve at ILD =1.8), and these effects 
are seen to be relatively small. The critical ILD of one 
radius is much lower than the Brown-Embury estimate, 
and it corresponds to a strain of 3.47% (f 1 ^ 3 — 0.15) for 
ILDo=1.8 at e = 10 9 /sec, corresponding to frame (c) of 
Fig. In the very final stages the ligament is drawn 
under biaxial stress, and the flow switches from radial 
material transport to tangential transport as the mech- 
anism switches from loop punching to drawing. At this 
point, the material is highly defective but it remains duc- 
tile. There is no abrupt fracture, as might be expected at 
larger length scales. Coalescence results from extended 
drawing and thinning of the ligament until rupture. 

Another measure of void interactions is whether the 
voids grow preferentially toward their neighbor. This ef- 
fect is quantified in Fig. 03 which shows the motion of the 
center of mass of the void surface (the void center) for the 
voids shown in Fig. Here ILDo=1.8 and e = 10 9 /sec. 
After the void growth starts, the center of void A ini- 
tially moves only slightly, but at about f 1 / 3 — 0.15 (ILD 
= 0.5 in Fig. [2j, it starts to move in the direction of 
the other void as the void growth becomes biased toward 
its neighbor. Just before coalescence the center of void 
A begins to move away from void B, as the growth is 
biased in the opposite direction. During this sequence, 
void B initially grows away from void A, then roughly in 
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FIG. 4: Growth of the voids until coalescence presented by 
void fraction for ILDo=0.5, 1.0, 1.2, 1.5, 1.8, and 4.6 diam- 
eters as well as in a single void case (in same box size) at 
i = 10 9 /sec. The asymptotic behavior (before finite size ef- 
fects) is exponential growth with exp(200e) as seen in the inset 
for the single void case. In the main figure the void fraction / 
has been scaled with the exponential. The circles point where 
the dynamical ILD's cross the line ILD= 0.5 in Fig- El 



unison with void A (/ 1//3 = 0.15) it begins to grow to- 
ward its neighbor, and before coalescence it too switches 
to growth away from the proximal void. This retrograde 
growth happens at the same point (after f 1 / 3 = 0.19) as 
the decrease of ILD begins to slow down in Fig. [21 [see also 
the snapshot in Fig. dd)]. The same phenomenon-first 
slow movement or repulsion from the void; then growth 
toward the nearby void after f 1 / 3 = 0.15 (ILD = 0.5); 
and finally retrograde growth-holds in the e = 10 8 /sec 
case, too, Fig. [3] As a reference the movement of the cen- 
ter of a single void in a box with e = 10 9 /sec projected to 
the same line is plotted. Comparing the single void case 
with the interacting voids at same strain-rates, one sees 
that the maximum distance the centers of the interacting 
voids have moved is three to five times larger than the 
nanoscale random walk of the single void center, so the 
movement of the void center is not just due to statistical 
fluctuations at the void surface. 

We have identified the onset of coalescence, but it is 
also interesting to examine the void growth following the 
onset of void interactions but prior to the actual coales- 
cence. How does this differ from the exponential growth 
of an isolated void 0, 0? In order to analyze the cor- 
related growth, we have factored out the non-interacting 
growth rate, exp(200e), from / in the plot in Fig. 21 The 
factor of 200 in the exponential is derived from the single 
void case, as indicated in the inset of the figure. The void 
growth data for ILDo=4.6 and 1.8 coincide with the sin- 
gle void curve. The void growth rate with smaller ILDo's 
reach their asymptotic growth rate earlier. In the figure 
we have drawn as circles the void size values, where the 
dynamic ILD's cross the line ILD=0.5 in Fig. O As can 
be seen from figure, there is no marked change in the void 



volume behavior when the voids start to interact. 

We have also performed a series of simulations of a 
fixed void in varying box size in order to find the coales- 
cence process of the void with its (six) periodic image (s), 
similar to the manner in which some continuum calcula- 
tions of coalescence have been done. The details will be 
reported elsewhere |20j, but it is worth noting, that the 
behavior is opposite to the results above in an important 
way. The smaller the box size, and hence the smaller the 
ILD, the later the void starts to grow. 

To summarize, interaction and coalescence of two voids 
in copper under tension have been simulated in multi- 
million-atom MD simulations. The effect of interac- 
tions between voids has been quantified by the increased 
reduction-rate of their separation and the movement of 
their centers. The interaction between the voids can be 
also measured in their shape changes |9| as will be done 
elsewhere 20]. The critical inter- void ligament distance 
has been found to be close to one void radius, indepen- 
dent of the strain-rate or the starting separation ILDo- 
The onset of coalescence occurs at the point that the plas- 
tic zones surrounding the voids first interact strongly. 

This work was performed under the auspices of the US 
Dept. of Energy by the Univ. of Cal., Lawrence Livermore 
National Laboratory, under contract no. W-7405-Eng-48. 
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